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Svjmmary 

The fundamentals of the Local Momentum Theory and a recent extension of 
this theory are presented in this report. The Local Momentum Theory is based 
on the instantaneous balance between the fluid momentum and the blade elemental 
lift at a local station in the rotor rotational plane. Therefore, the theory 
has the capability of e.aluating time-wise variations of air-loading and 
induced velocity distri.outions along a helicopter blade span. Unlike a complex 
vortex theory, this theory was developed to analyze the instantaneous induced 
velocity distribution effectively. This presentation of the fund;imentals I'f 
the Local Momentum Theory makes clear the boundaries of this theory and assists 
to make a computer prograan using this thei»ry. A new concept introduced Into 
the theory in this report is the effev?t of the rotor wake contraction in 
hovering flight. A comparison of this extended Local Momentum Theory with a 
prescribed wake vortex theory is also presented. The results indicate that 
the extended Local Momentum Theory has the capability of achieving a level of 
accuracy similar to that of the prescribed wake vortex theory over wide range 
variations of rotor geometrical parameters. It is also slnnm that the analyti- 
cal results obtained using either theory are in reasonable agreement with 
experimental data. 

Introduction 

This paper has two purposes; (1) to give a more (’■’tailed explanation vf 
the Local Momentum Theory than has been given previous papers^ and 
(2) to present recent improvements of the theory. 
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Thera ere two theorlei comonly used to eetlmete the Induced velocity 
distribution and the airloading distribution of a helicopter rotor, vortex 
theory and momentum theory. Vortex theory has the capability of analysing 
the induced flow field exactly, but the numerical integration requires 
considerable computer time. In addition, difficulties exist in the analysis 
of complex problems such as ptedicting rotor airloading under maneuvering 
conditions, even if a rigid wake system is assumed. 

In contrast, tlie ordinary momentum theory is simpler, being based on 
the assumption that tho induced velocity is uniform over the entire rotor 
disk'*. This, however, leads to unreliable airloading calculations. To 
overcome this shortcoming, various changes have been suggested. First, a 
rotor disk was assumed to be divided into annular ring elements'’, and 
momental balance was eonsidered on these elements. Tlien, the concept of 
the annular ring elements was extended to the pie-shaped elements*' in order 
to analyre the axlmuth-wise asymmetric flow field. Using these elements, 
a few analyses were presented tor a low frequency respoiase of a liellcopter 
rotor In addition to the nonuniformity in space, the time-wise 
variation of the induced velocity has been considered within the frame 
of the momentum theory. The eoncept of apparent mass was proposed for an 
entire rotor disk, the value ol which was determined by an experimental 
study of the thrust variation following a rapid collective pitch change**''’. 
The concept was then extended to apply to the pie-shaped elements, and im- 
provement of the predicting capability was achieved''”"’ for low frequency 
response of a rotor. As a further extension, the Local Momentum Tiieory 
was proposed by Ur. Azuma and this author'*-’ . 


In the Local >k>nentua TiTeory, a rotor blade is represented by a sories 
of n overlapping wings of decreasing wlng*-span. EUich overlapping wing has 
an elliptical circulation along its wing-span, producing, therefore, a uni- 
form downwash velocity. The actual airloading and induced velocity distri- 
bution of a real rotor blade is represented by the summation of the lift 
and induced velocity of this aeries of n overlapping wings. Momentum theory 
is used to obtain the relation between the lift distribution and the induced 
velocity distribution for each overlapping wingg, Consequently, in this theory, 
the momentum generated by a blade itself is counted. In contrast, in the 
previous momentum theory, the momentum generated by an entire or partial 
rotor disk is considered. Because of this basic difference, the theory was 
named Local Momentum Theory, ami it has the capability of predicting the 
higher frequency loading ol a helicopter rotor and also a more precise air- 
loading distribution along a blade span than Imve the previous momentum 
theories. The present theory, however, requires the use of a computer; tliere- 
fore, the theory cannot give an analytical solution wliich is useful in ob- 
taining a general view of a phenomenon. 

By comparison with vortex theory, the Local Momentum Tlieory can achieve 
iin equal level of accuracy, using a much smaller amount of computer time, 
llie reasons are as follows: the theory does not require the iterative pro- 

cedure nor the calculation of the inverse matrix; and numerical integratlen 
is avoided by using the attenuation coefficient, which expresses the time- 
wise variation of the induced velocity. The present theory, however, requires 
more complicated analytical procedure than that which is required by vortex 
theory . 
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Thtt thttory prasantsd In this paper la applicable not only to steady 
state conditiona, but also to various kinds of low frequency rotor responses, 
where trailing vortices play the primary role in the creation of the un- 
steady phenomena. The present theory is, however, based on steady airfoil 
characteristics: therefore, some extension of the theory is necessary when 
applying it to a problem of very high frequency rotor response, such us 
flutter, where shed vortices are the primary cause of the unsteady phenomena. 

An example of the extension of the present theory is shown in Ref. 21, where 
the unsteady airfoil characteristics, such as Wagner function, Ktlssner func- 
tion or the apparent mass for a blade, have been introduced in order to 
study the rotor gust response. 

Tlie Local Momentum Tlieory was recently extended to include the effect 
oi the rotor wake contraction in hovering flight. The theory previously 
assumed that the wake of a hovering rotor did not move in the radial direction; 
therefore, the theory gave alm''3t the same result as the rigid wake vortex 
theory'**, it has been pointed out, however, that the wake deformation 
has strong effects on the airloading and induced velocity distribution in 

hovering flight^*^*' ^ This extension of the theory to include the wake 

deformation effect was necessary in order to analyze airloading distributions 
measured by a laser velocimeter^ . In the new contraction model to be 
described in this paper, fluid elements surrounding a rotor can be moved as 
observed in experiments^ ^ In addition, the model to calculate the 
attenuation coefficient was improved in order to Include the effect of the 
blade tip vortex. The spanwise lift distributions obtained by this extended 
Local Momentum Theory are compared with experimental results with those of a 
prescribed wake vortex theory^ 
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Notation 


(Aij) 

AE 


a 

B 
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C* 
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1 
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K 
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k, 
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»m) 
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% 


M. 


aatrlx 

aspect ratio of a wing 

lift curve slope (rad”^) 

speed of sound 

constant or tip loss factor 

wing span or number of blades 

constavc or attenuation coefficient 

equivalent attenuation coefficient given by equation (57) 
thrust coefficient ■ 
wing chord 

function given by equations (15) and (lb) 
function given by equation (A, 1-9) 
function given by equation (A. 1-5) 
function given by equation (A. 3-7) 

rR i 

moment of inertia of a blade about flapping hinge - I (r-ITgJ drrv 

running index or inclination angle of tip path plane (positive 
forward tilt) 

coordinate of a blade element 
running index 

parameters showing tip vortex positions 
spring stiffness at flapping hinge 
1 if t 

airloading 

coordinate of a station on rotor rotational plane 
Mach number 

apparent mass of impervious disk given by equation (65) 

rR 

mass moment of a blade about flapping hinge - V f(^r- 
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» 

a 


n 

P 

q 

R 

r 




T 

t 

U 

V 


UMB of air 

aaat of air aaaociatad with local momantuw givan by aquations 
(21). (26) and (47) 

matbar of spanwisa partition 

rolling angular valocity of a rotor (positiva advancing-slda 
down) 

pitching angular valocity of a rotor (positive nosa up) 

rotor radius 

radial pcsiclon 

flapping hinge offset 

thrust 

time 

horizontal inflow velocity of a blade element givan by 
equation (29) 

frjrward velocity 


Vi.c 


Vn 




(X.Y.Z) 






y 

z 




mean horizontal inflow velocity of the i-th elliptical wing 
given by equation (33) 

normal component of inflow velocity (positive downward) 
induced velocity (positive downward) 
coordinate of fixed wing shown in Fig. 1 
nondimenaional radial position ■ 

nondimensional radial position given by equation (44) 

nondimensional radial position of tip vortex given by equation 
( 66 ) 


nondimensional flapping hinge offset « Vh 
spanwise position of fixed wing 

distance from rotor rotational plane (positive downward) 

nondimensional axial position of tip vortex given by equation 
(67) 
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Subscripts 

b 

cy 

down 

elp 


fUppitig anil* (potitiv* blad* up) - p§ + COStj^ + 

prcconlng anfl* 

circulation 

lock nundiar » 

graall incr«m«nt 

function given by equation* (50) and (53) 
nondinensional spanwis* position 

nondimensional spanwise position of fixed wing given by 
equation (6) 

blade pitch angle (positive leading-edge up) • 

•t ©.* s.'nV'-f •• • 

blade twist rate 

collective pitch angle at x «* 0.75 
running inde:: or inflow ratio 
advance ratio ■ 

nondimensional spanwise position of each elliptical wing 

air density 

summation 

solidity ^%.R 
inflow angle or wake age 
skewed angle 
azimuth angle 
rotor rotational speed 

number of blades 
cylinder 
downwash 
elliptical ■ ing 
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i 


J 


k 

P 

root 

tip 

up 

r 

X. 

0 


1c,1b 


0.75 


Suporicript 

(•) 


( ), 




■pamrif* partition, tha qu^atitjr of 1-th alllptlcal wins* 
or tha quwitity of tip vortax Mo. 1 

asinuth-iiriaa or tlna-wlsa partition, apamrlBo poaltlon, 
or tha qvumtlty of J-th alUptlcal trlng 

tha quwitlty of k-th blada 

rapraaantatlva point 

blada root 

blada tip 

upwaah 

circulation 

ailmuth-wlsa or tlna-wlaa partition 
uniform or initial valua 

first harmonic contants of Fouriar corina and sina sariaa 
tha quantity at thraa-quartar radial position 

d ()/dt 

tha quantity at a station (i, m) at tima t ■ J 


FundaMtntals of Ch« Local HoaMntun Thaory 



Fig. ahowa tha air flow flald aurro'muin(> an alllptical wing flying 
with conatant forward apeed V, A wall-known ^c^uU fro« tha lifting Una 
thaory^*^ ia that tha inducad velocity ia uniform inaide tha wing apan and 
the uniform induced velocity, 1* developed to 2V» far downatraan. Outaide 
the wing apan, there are upwaah regions, but the» are concentra;ad near the 
wing tips. Conaidering the taowentura change, the total lift of wing L ia 
given: 



u; 

where , ra « /* V 7C 

(2) 

Hie lift distribution, v , and induced 

velocity distribution, i/" , are: 


(3) 

3 iTo 

iT , 

f fcke K/i*»| spftH. (4.) 

[ 31/0 ( 1 - ) 

; l»fl out^side tU CS) 

where 



C6) 


Wlien a wing haa an arbitrary planfona or lift distribution, we may 
consider that the wing is composed of a series of n hypothetical elliptical 
wings. The elliptical wings may be arranged arbitrarily, for example. Fig. 

2^ shows a symmetrical arrangement and a one-sided arrangement. In each 
elliptical wing, tae momentum balance is given by equations (1) and (2). 

Using a ;?uffix to designate quantities of the i-th elliptical wing, equations 
(1) and (2) are rewritten as follows: 
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C7) 


where, 

TTiu* FVn (8) 

Ttie present theory was developed with the use of a computer; therefore, 
expreHaions with suffixes are adopted in this report for convenience in 
progr.iiwnlng. 

The lift distribution, li'fj), on the most left hand wing in Fig. 2' is 
expressed by the summation of elliptical lift, ti(0, on the coraponf'nt wings. 
As the lift distribution of each elliptical wing, ti(4), can be obtained by 
using a relationship similar no equation (3), the following equation is ob- 
tained ; 

hH) - s ^ (\ - g* W) 

where 

I » Ci “* 

and yoi shows the origin of the t. axis, that is, the mid-span of each 
elliptical wing. Substituting y - (-^ >7 from equation (6) into (10), 
equation (10) becomes 
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w 


whara 

%!,- (J2) 

In tha abova aquations, C is tha non-dlKanslonal coordinate fixed to 

each of the elliptical winga and J is the number of elliptical wings which 
include the point, ^ , within their span; therefore, only 
considered in equation (9). 

Substituting equations (7) and (8) into (9), the lift distribution, 
t can be exprewsed directly by the dowttwash of each elliptical wing 
aL/T, then 

03 ) 

llie induced velocity distribution of the most left hand wing in Pig» 2*, 

^nnh is also given by the summation of downwashes and upwashes of the component 
elliptical wings. As the induced velocity distribution of each elliptical 
wing is given by equations (4) and (5), the following can be obtained: 

+ Zf‘Vrj^(S) 04-) 

where 

= fufCS)-l- i 


j* 

f 

m 

f 

m- 
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In th« above aquatlona, ^down *xprasaas the uniform downwash Inside 
the wing span of each elliptical wing and f^p expresses the i;q>wash outside 
the wing span of each elliptical wing. According to blade element theory, 
the lift distribution of a wing is: 

ov 

Substituting equation (14) into (17). equation (17) becomes 

08 ) 

Now. two expressions of lift dlstrlbutlcn. that is. equations (13) and (18). 
have been obtained using Tnomentum theory and blade element theory respectively. 
In both equations, the only unknown variable is the downwash generated by 
each elliptical wing, i, 


Equations to solve for from equations (13) and (IB) are derived as follows: 

The lift per unit span. Jtj. acting on a local segment spanned by 
) (cross hatched area in Fig. 3) is given 

‘ i-v 

where, and are the tip positions of the J-th elliptical wing, and the 
(j -f l)-th elliptical wing respectively (Fig. 3). and they are non-dlmsn- 
sionalized by using equation (6). Substituting equetlons (13) and (18) into 
(19) glvss two expressions for ij. Then, equating these expressions and 
i.ilng equation (11), the following relation is obtainadi 
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M - ) ‘*Vcii - %i) 

Z ^i' (^^) 

15 1 

where 

and Yoi,is given by equation (12). Tlie above equation indicates that the lift 
due to blade element theory at a local segment spanned by (- ) is 

balanced with the local momentum change. 

Vflien the wing chord C(y) and the geometrical attack angle ) are 

given in the simple forms, equation (20) is analytically Integrated, and 
n X n simultaneous equations to solve for ( i,^l '^Tl.) are obtained 

by changing j value from 1 to n in equation (20).* 


(Ai-j) 

fAir, ^ 
■ i 

=. 

1 ' 


lAir„./ ’ 

Im 


VJlien equation (20) cannot be analytically Integrated, C(^) and are 

considered uniform inside the local segment spanned by (— 

and the values at the middle point of the local segment may be adopted: 
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C2J) 

Of 



In this simpliilcatlon, equation (20) can be integrated analytically, 
but a larger number of spanwise partitions is usually necessary to achieve 
a similar level of accuracy. 

Instead of the Integral method as stated in equation (20), there is 
another way to obtain nK n simultaneous equations, l.e., equation (22) to 
solve for aV[,( 71-) - Selecting a rapresentatlve point, Vrt , Inside 

the local segment spanned by ( ” ^^ ; -^hi ) and substituting Vn 
into equation (13) and equation (18) respectively, the following relations 
are obtained; 

i>t 

where 
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(26) 


» fX'Kil-iii 

ec7*'7») 


Although equations (25) anJ (26) have simpler form than equations (20), (21), 
(23) and (24), the author prefers the integral method (equations (20), (21), 
(23) and (24)) to the other (equations (25) and (26)). His previous cal- 
v'ulatlons for a rectangular wing showed that a smaller number of spanwlse 
partitions was required in the integral method than in the representative 
point method to achieve a similar level of accuracy. 

As shown in Fig. l’, the upwash regions outside the elliptical wing are 
concentrated near the wing tips; therefore, let us assume the upwash velocity 
induced by Oiich elliptical wing can be neglected, that is, fup (/. ) » 0 in 
equation (16). Substituting equation (15) into (14) and using tlie assumption 
fup V.) • 0, equation (14) becomes 

i/c/)= (27) 

L»l 


Using equations (23), (24), and (27), equation (20) becomes 



(28) 
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where lo given In equation (21). In the above equation, the only unknown 
variable ia C* instead of Al/t, ( (U» «* in equation (20); 

therefore, changing the value of J from 1 to n in equation (28), the 
solution of from 1 to n can be obtained without calculating the inverse 
matrix or any iterations. In other words, the coefficient matrix in equation 
(22), (Ac^l-)* becomes triangular by neglecting the upwash outside each 
elliptical wing. The sample calculation of this procedure and recurrence 
equation to solve are presented in appendix A-1. After the value of 

has been obtained in the above procedure, the induced velocity distribution 
and the lift distribution can be calculated by using equations (27) and (17) 
respect ively. 

An example of the application of the present method for calculating 
lift and induced velocity distributions for a rectangular wing wi^h aspect 
ratio (AR) » 6 is shovm in Fig. 4^. The number of spanwise partitions is 
n - 50. Tlie solid line and the broken line show the results when using the 
Local Momentum Theory with a one-sided arrangement and with a symmetrical 
arrangement respectively. The chain line shows the results when using 
lifting line theory''^ . The discrepancy in the results between the present 
theory and lifting line theory can be eliminated completely by introducing 
the upwashes that have been neglected in the present calculation. The error 
due to the neglect of the upwash is less in the middle of the wing span than 
in the other parts for the symmetrical arrangement. On the other hand, for 
the one-sided arrangement, the error is less in the tip of the wing than 


16 


in the other parte. Therefore, Fig. 4^ ehowe that the error due to the 
neglect of the upwash is lee* in that area of each elliptical wing which 
ia conanon to every other elliptical wing. Also, it ia shown that the error 
ia not very much even in the worst case, because the upwash region is limited 
to an area near the tip of each ellipticai wing as shown in Fig. l'. Wlien 
a more preoise estimation is needed in using the Local Momentum Tlieory, 
the upwash generated by an elliptical wing can be calculated as an iterative 
procedure. Usually, however, this procedure is not necessary when applying 
the present theory to the rotary wing. 

Wlien trying to applv the present theory to a rotary wing, the following 
dtflerenees between the lined wing and the rotary wing should be considered; 
U) the lorward speevl of a rotor blade is not uniform along the blade span 
because the rotor blade Is rotating, .ind (.2) the following blade may pass 
through the local station on tlie rotor rotational plane where the induced 
velocity generated by the preceding blade still remains. Therefore, the 
time-wise vari.ition ot the induced velocity generated by a preceding blade, 
should be estimated when analyzing the airloading of the following blade. 

I'o account tor these difterences, the theory is extended as follows. As 
shown in Fig. ■)', a rotor blade is assumed to be operated in a sheared I low 
in the rotor rotational plane. The forward velocity of the sheared flow is 
given along the blade span as follows: 

L)(jt) = v Siny + 


17 


wh«r« V 


is th« forward spaad of tha rotor, and \j/ la tha aalmuth angla of 
tha rotor blada. 

It la assuraad that tha rotor blada la compoaad of hypothatlcal allipti- 
cal wln^a almllar to the caae of tha fixed wing. (For convanlance, an el- 
liptical wing meana a wing which haa elliptical circulation along tha apan, 
although auch a wing does not hava elliptical planform nor elliptical lift 
distribution because of its operation in sheared flow.) By this approxi- 
mation of the flow, the trailing vortices shed from each elliptical wing 
may be considered to be straight, perpendicular to the wing span and to extend 
to infinity; therefore, referring to the Biot-Savart relationship, the 
induced velocity distribution along the span of each elliptical wing is 
the same as in the caae of the fixed wing, i.e. uniform within the blade 
apan, and it is given by equations (4) and (5). 

If each elliptical wing is one-sidedly arranged as shown in Fig. 5^, 
the lift near the blade tip can be estimated more correctly than the lift 
in other parts of the blade by using the Local Momentum Theory, as discussed 
above. As the greater part of lift acts on the blade elements near the blade 
tip, this arrangement is especially suitable for the calculation of a heli- 
copter rotor blade, and it gives a better estimation of airloading. 

In each elliptical wing spanned by ( > R ), as shown on the right- 

hand side of Fig. 5^, the following momentum balance can be obtained (Refer 
to Appendix A-2) : 


u = 

(3o; 

•«-£.=• r u \/i,c 

(30 
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where 


(JS) 


bi,» R(l- Jii.) 

Vi.,c“ + XlR ( 1+ 5^0/2 


In the above equations, bj means the span length of each elllptioal 
wing and Vj expresses the mean speed of the sheared flow along the 
span of the elliptical wing (Refer to equation (29)). Referring to 
Appendix A-2, the lift distribution on each elliptical wing, « 

is given as follows: 


where U (X). b{ and Vj are given by equations (29), ()2) and (ii), 
respectively. ■ is the non-dimensional axis fixed to each elliptical wing 
and it is related to the non-dimensional rotor radius, i.e. x, as shown 
in the following equations, because a one-sided arrangement is adopted 
here (Refer to Fig, b*). 




( 3 « 


Similar to the case of the previous fixed wing, the lift distribution and 
the induced velocity distribution of the rotor blade shown on the most left- 
nand side in Fig. 5^ are given by the summation of those of the component 
elliptical wings. These concepts are schematically shown in Fig. 6^. 
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Using sqiuitlon (34) and naglactlng upwash outslda aach alllptlcal vrlng, tha 
following aquations ara obtalnad: 

d . . 

- (2r v'fcc l-S* 06 ) 

I/U )- OD 

where 1^ means the number of alllptlcal wings which Include the point, x, 
within their span. 

In this subsection, for convenience in numerical calculations, many 
variables are expressed by using subscripts. Generally, the first subscript 
of any quantity, 1, indicates the i-th radial segment of the blade, the 
second subscript, j, indicates the time or azimuth-wise location of the 
blade segment, and the third subscript, k, indicates the k-th blade of a 
b-bladed rotor. When using this subscript form of expression, variables 
may be considered as uniform inside a given blade segment. For example, 
the spanwise induced velocity distribution, V(x) in equation (37), is 
also the azimuth-wise variable and may have a different value for each 
blade; therefore, V(x) is expressed as *^hat is, the induced 

velocity generated by the (i, j, k) blade segment. However, for the sake 
of simplicity, more abbreviated expressions are sometimes used. In 
equation (37), for example, Al/L is used instead of , although 

AV is the spanwise and azimuthwlse variable and is also a function 
of each blade. 
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Th« other expression of the spsiwlse lift distribution nay be derived 
by using blade elesMnt theory as follows: 

where (p is the inflow angle. In the case of the rotary wing, the induced 
velocity generated by the preceding blade nay affect the inflow angle, How~ 
ever, for simplicity, the equation to obtain the lift distribution is dis- 
cussed first, and the mathematical model to estimate the Inflow angle 
will be presented later. 


From the above discussion of the "integral method" in the case of the 
fixed wing, we may now derive simultaneous equations to solve ^l/^from 
equations (36), (37) and (38). Here, the integral area is chosen as 

) t because and^j^^fare the tip positions of the 1^-th elliptical 
wing and elliptical wing respectively (Refer Fig. 3). Equating 

the lift per unit span.^j', acting on a local segment spanned by ), 

the following equations can be obtained; 

c j a-c' - <P‘-' ) - x^) 


where 



isfl) 



XlRXc',t 

m 

CiS « TCii,c) 

(41) 

QC « 0(X»Xi,^c) 

(« 

OiAi » (X' (X* Xt/.c) 

(43) 

Xi/,c-CXi'+ 

(44) 


(*f) 
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(’ 41 .^ 








In th« abov« •quation*, bi, 5 «nd iTij^k «r« glvan by aquation# (32), 

(33), (35) and (37) raapactlvaly. Aa dlacuaaad In tha casa of tha flxad wing, 
equation (39) can ba analytically integrated whan U(x), C(x), 6(x) and a(x) 
are given in tha sinqtle form. However, for convenience, expreaaiona with 
aubacripta are preaented here, and U(x), C(x),8(x) and a(x) are repreaented 
by valuea at the middle point, Xc^6 . of the local aegment (Refer to equationa 
(40), (41), (42), (43) and (44)). For convenience in numerical calculation, 
the azimuth angle l/i^iB expreaaed aa follow#: 

t 




m 


where ^ and are the initial azimuth angle of the k-th blade and the 
azlmuth-wiae increment of the numerical calculation, respectively. 

Similar to the case of the fixed wing, the combination of equation (37) 
and equations (39) through (48) gives the recurrence equation to solve for 
if is given as discussed below. The complete procedure and the 

recurrence equation are presented iu the ^pendix A-3. 

Now let us discuss the inflow angle, , In equation (45). The flow 
velocity normal to the ro- rotational plane is generally composed of the 
induced velocity and the normal component of the rotor forward speed Vy . 
Furthermore, in the case of the rotary wing, the induced velocity may be 
generated by the preceding blade In addition to the blade under consideration. 
For example, in the case of the forward flight of the rotary wing. Fig. 7^ 
shows that the air flow at a local station (i, m) on the rotor rotational 
plane is influenced directly by two blade passages, i.e., by the i-th blade 
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of the k-th bl«4« «t Cl»t j or (1, j, k) blado sagMtiC mnd hy ch« 
i-th blade elaoMnt of the k^th blade at time J' or (!', J', k') blade seg- 
Ment. Conaequentl/i the nomal component of the total flow velocity at a 
local station (t, n) on the rotor rotational plane is given by the sumnatlon 
of the following! 1) Vj|» the normal component of the rotor forward speed 
(given by equation (46)), 2) the Induced velocity generated by the 

blade under consideration, 3) , the entire remaining Induced velocity 

generated by blades that have previously passed through the local station 
(« , m). Among the throe kinds of flow velocity, only the induced velocity, 

! k» generated by the blade under consideration is related to the lift 
when momentum balance is considered (Refer to equations (36) and (37'^). 
Dividing the normal coi^onent of the total flow velocity by the level 
forward velocity of the blade element, equation (45) can be obtained. 

Here, let uh discuss the time-wise variation of the induced velocity 
on the rotor rotational plane. It is shown in Fig. 8* that the rotor rota- 
tional plane is assumed to be divided into small elements called stations. 

The position of each station is given by the coordinate, (i, ra) . A blade 
element is shown to proceed on the rotor rotational plane within a small 
time 's.lerval for the case of hovering or steady vertical climb. Tlie 
blade element is located at a station (e', ra') at time t - J-1 with forward 
speed of it moves to stations (i, m) and (£", m") at time t - j 

and t " j+1 respectively. At time t «■ j-1 shown in Fig. 8(a)^ the normal 
velocity at the station (£', m') is given by the summation of the three 
conqionents l.e., Vh> and ^ ® small time Interval has 
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puiad, th* disturbed *lr at ataticm (f, «') haa ton* dowiward and 
th# fiald In tha rotor rotational plana la partially fillad with fraah air. 
Tha inducad valoclty at station (£' , ■*) haa baan changad at tlaa t • j, 
therafora, whila tha normal componant of tha forward valoclty, Vj|» la still 
constant (Rafar to Pig. 8(bJ). By introducing tha tlna-wlsa changing rata 
of tha 5nducad valoclty, tha raaaining inducad valoclty, V^/ , at station 
(■', ra*) at tima t ■ j can ba ralatud to tha pravious inducad valocity i.a. 
and ,, as follows j 

Vt H If ) 

Tliis time-wise changing rate of the induced velocity, i* callad "an 

attenuation coefficient." Generally, each station has its own attenuation 
coefficient, and the attenuation coefficient is a function of time. 

Because the blade element induces the velocity only at the station 

where the blade element is just passing through, equation (49) can be 
expressed in a more general form as follows: 

* 1,1,^ i- «• 

where ^ 4^18 one if any blade element e;.i8tB at station (1, m) at t ■ j-1 
and otherwise it is zero. Thus, if the attenuation coefficient is 

known, a combination of equations (39-48) and (50) gives the solution for 
^]/l with a given initial condition and a specified blade pitch input. 

After the value is determined, the induced velocity distribution can 

be obtained by using equation (37). The lift distribution is given by using 
equations (37), (38) and (45). Methods to obtain the attenuation coefficient 
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wiii b« discussed in ths next subsection 


In the esse of the irotery uing, each assuaed elliptical wing is one- 
sidedly arranged; therefore, when suaadng up induced velocities of each 
elliptical wlngf^ upwashes are concentrated outboard of the blade tip as 
shown in Pig. Although each upwash outside an elliptical wing is negli- 
gible, the suaraatlon of upwashes outboard of the preceding blade Lip is 
sn«M times strong enough to give an effect un the lift distribution of the 
following blade (Refer to Fig. 9*). Consequently, it is necessary to take 
this into account when estimating the lift distribution of t..e following 
blade operating outside the tip vortex of any preceding blade. Remembering 
the calculation procedure to obtain the value of the induced velocity of each 
elliptical wing, i.e. , It is necessary to omit the upwash only inboard 

of the blade tip (x 1 . This is because in the one-sided arrangement, 
the angle of attack at any spanwise station of the blade Is irdependent of 
the tip-side upwash, unlike the inboard upwash. Tlierefore, after the value 
of the induced velocity, Al/c , Is determined by omitting the upwash inboard 
of the blade tip (x 1), it is a slnq^le matter to estimate the upwash 
outboard of the blade tip (x 1) as follows; 

VU>n = £At4(,- (J,; 

Substituting equation (35) into (51) and using the expression with 
the suffix, equation (51) becomes 

° ) (s2) 
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wh«r« 1 It larger than n because it located outboard of the Made tip. 
Thus, wh**n the blade under consideration it operating outside the tip vortOA 
of any preceding blades, (.nuation (52) should be added to the calculation 
procedure stated before. Furthermore, in equation (50), the Induced velocity 
■'T.t. at station (t, m) should be given by the summation of all velocities 
Induced, not only those Induced by the real blade elements (x |» 1, or 1 ^n) 
passing through the station, but also those Induced by the hypothetical blade 
elements (x > 1, or 1 ^ n) extending through the tip-side upwash region. In 
addition to this, similar to the downwash, the time-wise variation of the 
upwash is expressed by using the attenuation coefficient. Consequently, in 
this case, equation (50) becomes 


where n* is the number of the spanwise partitions of the real and hypothetical 
blade. considered to be one if any real or hypothetical blade 

element exists at station m) at time t •• j-1; otherwise it should be 
considered to be zero. 

When the pitching or rolling motion of a rotor hub or the blade flapping 
motion is considered, only the normal velocity component needs to be 
modified ao follows: 

VV j * V a>i ^ 
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It should bs rsAllsod that tha normal componmt Is now tha function of 
tlma or tha axlmuth angla in tha abova aquation. Similar to tha modifica- 
tion shown in aquation (54), more complex hub mutlon or tha blade deformation 
can also be introduced into the present chaory without any essential change 
of the inflow calculation procedure. 

The only difference between the present theory and other numerical 
calculations comes from the procedure to obtain the lift or induced velocity 
distribution from the given blade operating condition. Therefore, the cal- 
culation procedure is very similar to the usual computer provr jm used for 
rotor aerodynamic analyses. One sample calculation procedure using the 
present theory is shown in Fig. 10. 


As the present theory is designed to trace a time-wise variation of 
the lift or of the induced velocity distribution, a suitable numerical 
nkHliod, such as the Rung«.*-Kutta method, should be used in order to solve the 
blade flapping equation. As an examiile, the sample blade flapping equation 
is: 



where is spring stiffness at the blade root and is the preconlng or 
neutral angle. Tlie shear force at the blade root can be obtained by the 
numerical integration of the lift distribution from the root to the tip 
along the blade span, and the rotor thrust is obtained by the summation ol 
the shear force at each blade root. 
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Becausa th« initial inducad valocity dittrlbutlon and tha initial blade 
flapping motion are unknown, a calculation to obtain tha initial trtead 
condition is required at the beginning of the computation. In tha case of 
the calculation of a steady state, this initial trimmad condition is the 
solution, and the computation is completed. This initial trimmed condition 
is usually obtained within eight revolutions of the rotor in hovering, and 
within five revolutions in forward flight 0.15). The initial value of 

the Induced velocity is assumed to be jsero in the present program. It is not 
effective in shortening the initial transient period to use the value of the 
Induced velocity instead of zero. This a because this initial 

transient period depends mainly on the number of blades, on the development 
1)1 the wake, and on the damping of the blade flapping motion. 

In the actual calculation using equation (5i), the following dis- 
cussion might be useful; as stated before, the axes [t, m] exist 
on the rotor rotational plane, and the origin is fixed to the Inertial 
sp.i. e in the case of a steady flight. Toe axes [i, J, klAr«flxed 
to the blade, and the origin exists on the hub. The small difference 
between tlie two reference frames due to the coning angle is Ignored. In 
the case of hovering or vertical climb, the axes [I, j, k] haveouly a 
rotating motion relative to the axes 1«, m); therefore, the pie-shaped 
element is recommended for the station (t , m) , for the convenience ol 
the conversion between the two reference frames. In contrast. In the 
case of the forward flight, the axes [i, j, kj have the transversal motion 
In addition to the rotating motion relative to the axes l^, ml; therefore, 
in this case the square element is preferred for the station (t , m). 
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Atfniuition Coefficignt 


As statad in tha pravlous sactlon, tha attanuatlon coafficiant is Intro- 
ducad into tha prasant thaory in ordar to axptass tha tima-wise variation of 
the Induced velocity at a local station on the rotor rotational plane. When 
it is necessary to estlmata the induced velocity distribution on another 
specified plane such as the tail rotor operating plane, the attenuation co- 
efficient should be obtained for that plane using a model which is different 
from the one presented here. The information obtained for the induced velocity 
distribution on the rotor rotational plane is usually enough, however, for 
an analysis of the rotor airloading. Consequently, discussion in this paper 
la limited to the use of the model in obtaining the attenuation coefficient 
on the rotor rotational plane. 

Tlie attenuation coefficient can be defined as the time-wise changing 
rate of induced velocities between any arbitrarily chosen time and that time 
when the blade element passes through the local station (e, m). Therefore, 
if the exact value of the attenuation coefficient is desired, the physical 
model used should be very complex, such as the distorted wake vortex theory. 
Such a complex model requires greater computer time, however; consequently, 
the balance between accuracy of calculation and complexity of the model is 
the most Important consideration in determining the model to be used in 
obtaining the attenuation coefficient. 

Let us assume steady flight of hovering or of vertical climb, and, 
remembering the calculation procedure stated in the previous section, let 
us limit the purpose of the present analysis to the estimation of the lift 
distribution acting on a rotor blade. It is necessary, then, to know the 
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induced velocity rameining at station (I, in) only if after the preceding 
blade has passed, the following blade is located just before the station. 

The attenuation coefficient at station (£, m) can be given, then, as the 
ratio of the Induced velocities of two times; one, when the preceding blade 
is passing through the station, and the ether, when the following blade is 
located just before the station. The period between the two moments is 
given by 

Atir® 

The period aU, in the above equation, can be obtained before the beginning 
of the calculation, and it is constant during the entire calculation pro- 
cedure; therefore, during the calculation, the attenuation coefficient is 
no longer a function of time but a function only of the position of the 
station. 

Furthermore, in the case of the simplest model, the attenuation 
coefficient is assumed to be uniform over the entire rotor disk. This 
simplest model was used in the earliest stage of the development of this 
theory to check the characteristics o£ the attenuation coefficient. The 
valte of the uniform attenuation coefficient is determined by comparing the 
theoretical thrust coefficient with the experimental results. Thrust co- 
efficients of a rotor were obtained theoretically by using the Local Mo- 
mentum Theory, assuming various values of the uniform attenuation coefficient 
selected between 0.5 and 0.95. The curve of the thrust coefficient versus 
the uniform attenuation coefficient for a given collective pitch is shown 


u a solid line in Fig. 11^. The thrust coefficient tends to become small 
when the value of the unifoxrm attenuation coefficient is increased, because 
more Induced velocity remains on the rotor rotational plane and the angle 
of attack of the rotor blade decreases i Eventually, the thrust coefficient 
becomes sero at the point C “ 1,0, where all induced velocity remains; there- 
fore, new Induced velocity to be balanced with the lift acting on a blade 
is added to the previous induced velocity until the angle of attack of the 
blade becomes zero over the entire disk. 


Tlie chain line In Fig. 11^ shows the thrust coefficient obtained by 
an experimental study^ using a model rotor. This thrust coefficient, of 
course, is independent of the attenuation coefficient. The intersection of 
theoretical (solid) and the experimental (chain) lines gives the value of 
the uniform attenuation coefficient for a given collective pitch. 


In Fig. 12'^, the intersections for various kinds of model rotors 
are*^*'^^"* plotted. As the present uniform attenuation coefficient is 
defined by the period given by equation (56), the attenuation coefficient 
depends on the number of blades, i.e., b. Consequently, the equivalent 
attenuation coefficient C* is used for the lateral axis in Fig. 12'' in 
order to unify various attenuation coefficients depending on the number 
of blade. As a three-bladed rotor is used as the standard in this figure 
(Fig. 12*’), the equivalent attenuation coefficient of a b-bladed rotor is 
defined by 



(sV 
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The next step In estimating the attenuation coefficient Is to develop 
a purely theoretical model. As shown In Fig. 13, the air mass hit by a 
preceding blade (Blade (A)) at a station (I, m) at time t - j-1 moves 
downward as time passes. VHien the following blade (Blade (B)) comes to 
the same station (1, m) at time t - J, the distance between the air mass 
and the station (l, m) Is shown as Z In Fig. 13. Here it is assumed that 
the air masses inside the rotor disk move downward with the same constant 
speed, in other words, the circular disk composed of the air masses moves 
downward. Such a flow field can be expressed by a semi-infinite vortex 
cylinder of uniform disk loading, and the circular disk moving downward 
corresponds to the upper end of the cylinder. Because the flow around the 
vortex cylinder is determined by using simple vortex theory^ the attenu- 
ation coefficient is given as the ratio of the induced velocity at the rotor 
rotational plane and that at the upper end of the vortex cylinder, that is, 
between the two hatched elements shown in Fig. 13. After the blade (B) 
passes through the station, the process starts again. In other words, it 
is as if the cylinder suddenly returns to the rotor rotational plane and 
again begins to move downward. Each station, then, has its own vortex 
cylinder which moves up and down at each blade passage. Consequently, the 
attenuation coefficient obtained from the semi-infinite vortex cylinder is 
a function of two variables. One is the position cf the station, i.e., 
coordinate (e, ra) , and the other is the length between the upper end of the 
vortex cylinder and the rotor rotational plane that is shown Z in Fig. 13. 

In the case of hovering, Z is given by 

z = vT. Atb isS) 
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whartt Atir Is obtainad from aquation (56) and^/o doimward moving 

valoclty of tha uppar and of tha vortax cyllndar; tharafora, 

Tha theoretical attenuation coaffiant obtainad from tha saml-lnf Inlta vortax 
cylinder In the case of hovering Is shown in Fig. 14^. 

The merits of this theoretical model for the case of hovering are as 
follows. First, the attenuation coefficent is a function of three variables, 
the coordinate (t, m) , the length Z given by equation (58) and the span- 
wise position Therefore, a general chart of the attenuation coefficient 

can be obtained (Fig. 14-^) before the beginning of the calculation, and the 
chart can be used for any single rotor configuration. In addition to this, 
the flow field around the serai-infinite vortex cylinder can be obtained with 
very short computer time because of the available analytical solution'^. 

Using equations (56), (58) and (59), the attenuation coefficient can 
be calculated for each model rotor from the given thrust coefficient and the 
specified position of the station. The result at 0.7 and Vr - 0.8 is 

compared with the previous experimental results in Fig. 12^. The reason- 
able agreement between the two results verifies the theoretical approach 
presented here. Hie result shows that the theoretical value at 9k- 0.75 
may be used as the uniform attenuation coefficient when the simplest cal- 
culation is desired. 
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Th« use of the theoretical model in hovering flight to obtain the 
attenuation coefficient can be extend«l easily to the case of forward flight. 
Similar to the hovering case, the serai-infinite vortex cylinder with uniform 
disk loading is used to give the attenuation coefficient. In the case of 
forward flight, however, the cylinder is assumed to be inclined, the angle 
of which is given by 

m 

Hie flow around such an inclined serai-infinite vortex cylinder is known from 
H imp If vortex theory’ ^. 

Furthermore, the attenuation coefficient in the case of forward flight 
Is asauraed to be uniform over the entire rotor disk. This assumption is 
based on the following considerations. From a moving axis fixed to a 
helicopter, a station on the rotor rotational plane appears to move 
rearward at the same speed as the forward velocity of the helicopter. 

Thus, every station has come from the upstream and has gone again to the 
outside of the rotor disk after a finite number of encounters with blade 
elements. Consequently, it is expected that the attenuation coefficient 
under the usual operating conditions of forward flight does not play as 
impo**tant a role as under hovering conditions, and therefore, the simpler 
model might be sufficient. 

In the case of forward flight, in addition to the above assumptions, 
it is also necessary to assume that the period during which two successive 
blades pass through the same station, is given again by equation (56). 
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Exactly Bpaaklng, th« pariod !• no longer glvan by iVMSh a slapla ^xpraaalon. 
Tha period la ahortar on the advancing aide (# H 90 °) mid longer on the 
retreating aide (♦Vs 270°) than the period given aa equation (56), becauae 
the rotor movea forward. Eowever, aa atated above, the alnpler model might 
be aufflclent In the caae of forward flight, and the rotational apeed of the 
blade la much greater than the forward apeed under usual operating condition; 
therefore, equation (56) glvea a aufflclent approximation of the period. 

Thus, the downward dlatance of tha movement of the vortex cylinder 
during the period la given by 

Z *■ CV +• ceSL/ • ^ +■ Vo ) 

where la the mean Induced velocity. Also, the rearward distance is 
given by 


Xcy ~ \/COS 


(a) 


Combining the analytical approach^^ to solve the inclined semi-infinite 
vortex cylinder with equations (56), (60), (61) and (62), the attenuation 
coefficient enn be obtained at any station. Remembering the discussion in 
the case of hovering, the value at Vr “ 0.75 and ip * 90° or 270° may be 
used as the uniform attenuaton coefficient in the case of forward flight. 
Thus, similar to the case of hovering, the general chart for the theoretical 
attenuation coefficient in forward flight can be obtained with very short 
computer time. The result is shown in Fig. 15^, and it can be used for any 
single rotor configuration. As will be seen in examples in the next section, 
this uniform attenuation coefficient based on the above assumptions gives a 
sufficiently good estimate of the variation of the Induced velocity in 
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forward flight. It is, howavtr, rtcoaMndad that tha attanuation coafficiant 
at low forward flight (u 4 0.10) ba a fimction of radius and asinuth angla^^, 
similar to tha casa of hovarlng. 

In tha casa of unstaady flight, tha naan inducad valocityt/^ is a 
function of tlma. Tharafora, the valua of \To is givan by tha avaragad in- 
duced velocity along tha blade span at each Instant. Substituting tha valua 
of ^To into equation (58) or (61), and using the general chart, the attenua- 
tion coefficient can be obtained at each instant. The introduction of tha 
time-wise variation of xTo into the present theory is nearly equivalent to 
time-wise change of the pitch of the helikal wake in the vortex theory. 

Furthermore, when the rotor hub has unsteady motion, equations (58) and 
(bl) should be modified. For example, when the rotor hub is rolling or 
pitching, equation (58) becomes 


Kquation (61) becomes 

Z= j + \/cosi >‘ - RX. + pSiV'l/*) 

Tlvus, the wake deformation and the interaction between a wake and blade 
airloading, which are caused by an unsteady hub motion or an unsteady control 
input, can be easily introduced into the present theory. 
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Application of th« Thaory 


Hovarlng Rotor 

Th« lift distribution along a blada span of a hovarlng rotor vws cal- 
culated by the theory given in the previous sectlont Ir the calculation* 
the nuaber of the blade sparwise partitions is n ■■ 20, and the asiauth-wise 
increment is , The lat ter may seem too large, but it is allowed 

because only the steady state with no blade-flapping motion is concernedt 

An example of the reault is shown in Fig. lb'. It is compared with the 
result obtained from vortex theory* and also co^ared with the experimental 
result'*' using a model rotor in the wind tunnel. As shown in the figure, 
the result of the present theory using non-uniform attenuation coefficients 
is very close to that of the vortex theory and is reasonably close to the 
experimental result. Kven using the constant attenuation coefficient, the 
present theory gives a good estimation of the blade spanwise lift distri- 
bution of a hovering rotor. Tlie computer time of the present theory is at 
most about 1/10 of that of the vortex theory. 

Advancing Rotor 

Wlien a rotor is operating in forward flight, the effect of the upwash 
velocity observed outside the preceding blade-tip must, as stated before, 
be included in determining the angle of attack of the following blade. 

An example of the analysis is shown in Figs. (17), (Id), (19), and (20). 
In the calculation, the number of blade spanwise partitions is n “ 20, and 
azimuth-wise increment is ^ 10 deg. Tlie rotor rotational plane is 

divided into a net of square elements of dimension R/80, the number of which 
is i'ra * lfa0*320. The lift coefficient is assumed to be zero in the reversed 
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flow rttginn. Th« rotor it ucuwid to h«v« articulated rigid bladaa and 
the airloading ii coupled with the blade flapping motion. The epanwlae and 
azlmuth'-wiae lift dlatrlbutlona are shown in Fig. 17^ and in Fig. Id^ 
respectively, being coeq;>ared with experimental results^^. For the sake 
of comparison, the result obtained by using a vortex theory^* is also 
shown In Pig. 18(a). In addition to the lift distribution, the radial 
distribution of time averaged induced velocity over one rotor revolution is 
shown in Fig. 19*.. The result of the present theory is compared with that 
of a simple vortex theory^ Furthermore, the diatributi<-n of the angle 
of attack obtained by using the present theory is shown in Fig. 20‘*’ . 

In this figure, the broken line indicates the trace of the intersection of 
the blade under consideration with the tip vortex of the preceding blade. 

As shown in Figs. 17* and 18*, the correlation between the present theory 
and the experiment is reasonable, even though the constant attenuation co- 
efficient is used for sinqilicity of calculation. The correlation near the 
intersections of tlie blade under consideration with the tip vortex of the 
preceding blade (Refer Fig. 20) is not, however, as good as the rest. Tl-iis 
is because the wake deformation is not considered in the present calculation. 
Since the lift distribution at the outer part of the blade is more sensitive 
to the tip vortex positions, the difference between the present theory and 
the experiment becomes greater at this part. 

In addition to the effect of the latest tip vortex, the effects of the 
older tip vortices on the angle of attack distribution are observed in the 
result using the present theory (Fig. 20^*^). These tip vortices cause the 
variation of the spanwise lift distribution (Fig. 17*), and the high frequency 
oscillations of the lift (Fig. 18*). The use of the constant attenuation 
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co«fflcA«nt t«nds to «n{>hMlM thtto •ffccta. Thii !• b«cau«« th« variable 
attanuatlon coafficlant near tha Inttraactloaa of the tip vortlcaa with the 
blade ia amaller than the conatant attenuation coefflcientt When uaing the 
variable attenuation coefficient, the effect of the lateat tip vortex alone 
la remarkable 

Aa ahown in Fig. 18(a), the lift obtained by the present theory la more 
quickly changed than that obtained by vortex theory. This la because tha 
present theory la baaed on the steady airfoil characteristics and Ignores the 
shed vortices. 

Rotor Response due to Rapid Increase of Collective Pitch 

A transient response of the rot^ir thrust and a related induced velocity 
variation following a rapid collective pitch change have been studied in 
experimental tests'’* Analyses also have been conducted by Introducing 
an apparent mass associated with the rotor disk into the classical momentum 
theory^* Such an analysis, however, requires the use of experimental 
results to obtain the value of the apparent mass. The approach using the 
vortex theory is orthodox but it requires much computer time, even though the 
rigid wake is assumed^’’. In contrast, it is very easy to apply the present 
theory to this problem and the result can be obtained with much less computer 
time. 


Fig. 21^ compares an analysis for the case of hovering using the present 
theory with one using the vortex theory^^. The aerodynamic thrust shown 
in the figure means the thrust which does not include the inertia force due 
to the blade flapping motion. 


39 


In th* calculation, tha maan Inducad valoclty , which la uaad In 
aquation (58) In order to calculate the moving dlatanre of the vortex cy- 
linder, la d<^termlned by two methodat by averaging the Induced velocity 
along the bicde apan at each Inatant, and by ualng tha aame value uaad In 
the calculation for the vortex theory. The aolld line ahowa the reaults 
ualng the former, and the chain line ahowa the reaulta ualng the latter. 

In eitaer way, the reaults of the present theory are coincident with those 
of the vortex theory «hown by the broken line. This indicates the insen- 
sitivity of the transient shape of t/s to the thrust coefficient. 

In the variation of the aerodynamic thrust, the first decrement following 
a peak results from the blade flapping motion which reduces the blade attack 
angle. The second decrement following the maximum thrust is due to the 
development of the wake which increases the induced velocity. 

A comparison is shown in Fig. 22^ of the results of the present theory 
and those of the classical momentum theory in which the apparent mass of the 
entire rotor disk has been introduced. The value of the apparent mass as- 
sociated wil i ' hi‘ rotor disk is assumed Mq or 2 Mq, where M is the apparent 
mass of the ^.npervious disk in the unsteady translation perpendicular to 
its plaii Mq is given by 

= (6S) 

Tlie blade is constrained at the flapping hinge in this calculation in 
order to corapare only the aerodynamic characteristics obtained by using two 
theories. The irregular curves of the thrust and the induced velocity ob- 
tained by using the present theory result from the mutual interference among 
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the induced velocities of all blades. In the calculation, the nuadier of 
blade spanwise partitions is n ■■ 20, and asirauth-wise increment •• 10 

deg. The present calculation supports Carpenter's experimental result^ ° 
tlukt the apparent mass of the rotor due to rapid collective pitch change is 
equal to that of the in^ervious disk in the normal translation. 

Rotor Response due to Sudden Increas e of Cyclic Pitch 

Fig. 23' shows the time history of the inclination of the tip path plane 
following the step cyclic pitch input into a hovering rotor. Tlie lines show 
the truces of the point the unit vector which stands at the hub center 
normal to the tip path plane. Tlie inclination of the unit vector is obtained 
from the individual blade flapping angle by using the multiblade coordinates* 
at each instant. As the sample rotor does not have flappin hinge offset nor 
a spring at the blade root, the coordinate of the steady state is (0, -1). 

Tlie solid line shows the results of the present theory. The broken line 
shows the result using the model of constant induced velocity .istribution 
over the entire rotor disk^'’. The double chain line indicates the result 
obtained by the theory^ based on the momentum balance in tlie pie-shaped 
area in the rotor disk. The results obtained by the introduction of the 
apparent mass of the rotor disk to the theory, in addition to the former 
tlieory itself*’, are shown by the chain lines. Tlie results indicated by 
the chain lines are reasonably close to those of the present theory for 
various values of Lock number. 

In the classical momentum theory, the time-wise variation of the 
induced velocity is approximated by the inertia term due to the apparent 
mass of the rotor. Therefore, the difference between the chain line and 


the double chain line shows the effect of the time-wise variation of the 
Induced velocity. Purthermore, the difference between the chain line and 
the broken )ine Indicates the effect of the space wise non-uniformity 
of the induced velocity distribution, in addition to the time-wise 
variation. As shown in the figure, the latter difference is small at the 
beginning of the motion, and the former difference becomes small as time 
passes. 

Tlie reasons are as follows. In the model shown by the broken line, the 
Induced velocity is constant; therefore, when the imbalance of the flapping 
moment la produced by the step cyclic pitch input, only the blade motion is 
changed In order to satisfy the new equilibrium. In the model shown by the 
double chain line, the induced velocity distribution can be changed ir addi- 
tion to the blade flapping motion. However, the induced velocity distribu- 
tion is given by the condition at the terminal steady atate^; therefore, 
tlie induced velocity has developed to the terminal value even at the very 
beginning of tlte transient motion. In the model shown by the chain line, 
due to the inertia term associated with the apparent mass, the induced 
velocity distribution is not changed as rapidly. Consequently, at tlie 
beginning of the transient motion, the Imbalance of the flapping moment 
is satisfied largely by the change of the blade flapping motion. Tlie 
induced velocity, then, develops to the steady state, in the case of 
the very low frequency response of the rotor, therefore, the effect of the 
apparent mass may b«; neglected, and the model shown by the double chain line 
can give a good approximation of the rotor inclination'^. 


# 
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Th« step cyclic pitch response of the advancing rotor is shown in 
Fig. 24^. The time histories of the thrust coefficient and the coning 
angle are indicated in Fig. 24 (a) The inclination of the tip path plane 
is shown in Fig. 24(b) The small high-frequency variation observed in 
the thrust coefficient results primarily from the non-uniformity of the 
Induced velocity distribution, caused primarily by the tip vortex of the 
preceding blade. Tlie effect of the inclination of the tip path plane on 
the thrust is shown as the low frequency variation in the figure. In the 
calculation, the number of the spanwise partitions is n - 20, and the 
azimuth-wise increment is ■ 10 deg. 

A few more applications of the Local Momentum Theory have been tried. 
For example, the tandem rotor in yawed flight and the coaxial rotor are 
analyzed in Ref. 31. The gust response of a single rotor is studied iu 
Ref. 21. 


Tlie new hovering model including the wake deformation 
Calculation model 

A recent extension of the Local Momentum Theory to include the wake 
deformation of steady hovering flight is presented in this section. Although, 
as discussed in the previous section, the results using the Local Momentum 
Tlieory showed good correlation with a few experimental cases, this extension 
of the theory was necessary in order to analyze the airloading distributions 
measured by a laser velociraeterff^*^ 

The wake deformation of a hovering rotor is separated into spanwise 
deformation and axial deformation. The spanwise wake deformation includes 
the contraction and rollup of the vortices. The axial wake deformation is 
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the axial deviation from the rigid helical wake poaition. In the Local 
Momentum Theory, the spanwlse wake deformation la equivalent to the apanwlae 
movement of a local atation on the rotor rotational plane. In order to 
expreaa the axial wake deformation, a new model to give a corresponding 
attenuation coefficient was developed. Because of the effect of the tip 
vortex on a spanwise lift distribution is dominant among all the wake com- 
ponents, the model of the present theory was developed primarily to trace 
the behavior of the tip vortex. Fig. 25(a) and Fig. 25(b) show the side 
view and top view of the tip vortex. Tlie tip vortex No. 1 in Fig. 25(a) 
is shed from the blade under consideration (blade (A)) at position at 
time t ■ j . The tip vortex No. 2 was shed from the preceding blade (blade 
(B)) at position Pj at time t ■ J-1, and is located at position P 2 at time 
t » j. The period between t * j-1 and t * j is the period during which two 
neighboring blades successively pass through the same position of the rotor 
rotational plane, that is, At-b given as equation (56). During the period 
Atb , the tip vortex No. 2 moves from Pj to P 2 * Following the motion of 
the tip vortex No. 2, a station on the rotor rotational plane moves from rj 
to r 2 during the period, where r^^ is a radial position of the tip vortex 
No. i. Similarly, each atation moves from ri_|to during the period 
following the motion of the tip vortex No. i. The positions of the tip 
vortices are given by a generalized wake model^^^* ^ 

Xt= K+^-C(-K+)e- (f/h) 


= K.^ 
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vrh«r« and ara the axial and radial non-dimensional psiaitions of the 
tip vortices. Using aquation (66), the radial position of the tip vortex 
No. i is given as 

Outside the rotor disk (r > R), each of the stations moves a distance equal 
to that between r^ and V 2 toward the blade root In order to avoid making a 
gap between neighboring stations. At the inner area of the rotor (r^RK^» 
the stations don't move, because no tip vortex moves into this area (Refer 
to equation (bb)). Tlierefore, the positions of the stations in the inner 
area are coincident with those of the rigid wake model. The parameters 
K 2 , K 3 and K 4 in equations (bb) , (b7) and (b 8 ) are determined by experi- 
mental data using flow visualization techniques. There are three available 
generalized wake models, one proposed by Landgrebe^^’, ,one by Kocurek and 
Tangular^'', and one by Kocurek, Berkowitz and Harris^-’. In this paper, 
the second model is most frequently used, because of the applicability over 
wide range variations of rotor geometry, and because of the convenience of 
comparison with other computer codes. 

As each station moves, a part of the following blade (the hatched area 
in Fig. 25(b)) operates in the upwash region generated by the preceding blade, 
even in the hovering case. Tlierefore, it is necessary to calculate the up- 
wash distribution outside the blade tip in addition to the downwash distri- 
bution along the blade. The upwash or downwash distributions are calculated 
at each station by using equation (37), equation (39) through (48) and 
equation (52). 
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The new model was also developed to take the wake deformation into 
account and to obtain the attenuation coefficient at the upwash region. 

Fig. 26 shows the new model which Is composed of a vortex ring and a seml- 
Inflnlte vortex cylinder. The vortex ring represents the tip vortex nearest 
to the rotor rotational plane, and the vortex cylinder represents the re- 
maining tip vortices. The vortex cylinder is made of the continuous vortex, 
the strength of which la equal to that of the axially averaged tip vortices^^. 
The strength of the vortex ring is equal to that of the tip vortex. Fig. 

26(a) shows that the preceding blade passes through a local station at time 
t » j-1. The vortex ring is localed on the rotor rotational plane, and the 
top of the vortex cylinder is located at the position of the second tip 
vortex. The distance between the top of the vortex cylinder and the rotor 
rotational plane is Z 2* ^^8* 26(b) shows that the following blade comes 

lust before the station at time t • ,j . During the period between time 
t “ .1-1 and t - 1, the vortex ring and the vortex cylinder move downward 
and radially contract. Using equation (67), the vertical position Zj is 
given by 


= RK,#>, 




CIO) 


where 4>j is given by equation (69). Tlie radius of the vortex ring is equal 
to that of the tip vortex nearest to the rotor rotational plane (Refer to 
Fig. 26(b)). For simplicity, the radius of the vortex cylinder is assumed 
to be equal to that of the vortex ring instead of the corresponding tip 
vortex radius, This leads, however, to an overestimation of the radius 
of the vortex cylinder because of the contraction of the tip vortices as 
shown in Fig. 25(a). The radii of the vortex ring and the vortex cylinder 
are given as equation (68). 
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The attenuation coefficient is defined by the changing rat« of induced 
velocities at a station (t, m) between time t - j and t - j-1. In the 
present model including the wake deformation, the station (t, m) moves from 
position rj_^ to position r^ on the rotor rotational plane during this 
period. Consequently, the attenuation coefficient of that station is given 
as 

ctl = 

where is the induced velocity at position r^_j^ on the rotor rotational 
plane at time t - j- 1 , and is the induced velocity at position r^ on 

the rotor rotational plane at time t ■ j. Because the station at the inner 
pari i>f the blade ( riRK4) does not move, the attenuation coefficient of 
that station is given as the changing rate of the Induced velocities at the 
same radial position between time t “ j and t •• j-1. The stations which 
are located between positions T 2 end R at time t » j were situated 
outside the rotor disk (r ' R) at time t “ j-1. Therefore, when determining 
the attenuation coefficients of these stations by using equation (71), the 
induced velocity is calculated outside the rotor disk. If the old 

model, which is composed of the vortex cylinder alone (Refer to Fig. 13), 

Is used, the induced velocity is zero outside the rotor disk, and it 

Is impossible to define the attenuation coefficients of these stations. 

Consequently, the discrete vortex was introduced into the present model. 

In the Landgrebe wake model and in the Kocurek and Tangular wake model, 
the parameters Kj^, K 2 , and K 3 in equations ( 66 ), (67), ( 68 ) and (70) are 
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funcclons of a Chruat coafflclant^^* In tha prasant mathod, thosa 
parametara ara calculatad at aach Lima atap by ualng tha value of a thruat 
coafflclant. Thera la an analytical aolutlon^^ for the Induced velocity 
around the vortex cylinder or the vortex ring. Baaed on th.'a solution, 
general tables of the Induced velocity are constructed for the vortex ring 
and for the vortex cylinder, before the start of the calculation. The 
attenuation coefficient is calculated at each time step from these general 
tables by using equation (71). The entire calculation Is completed when the 
thrust coefficient reaches a steady state. The calculation procedure of the 
present theory is shown by Fig. 27. 

Calcula tion Results 

A comprehensive study has been conducted to check this extended Local 
Momentum Theory. Tlie spanwise lift distributions for five different rotors 
were calculated by this theory. The results were compared with experimental 
data and with those obtained by using other analyses. Table 1 shows the 
geometrical characteristics and the operating conditions of sample rotors 
used in this study (see also Refs. 26, 34, 36-38). Three of them are model 
rotors, and two are full scale rotors. 

Two computer codes of the prescribed wake vortex theory were used in 
this study, the MI lifting surface program'’^’, and the UTRC lifting line 
program^ In the present theory and the UTRC code, the lift curve slope 
of airfoil sections is given as 




a2) 

do ~ 


(7i) 
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M !• th« Mach niuobar of tha airfoil faction and it ia givan aa 

M - ilR Xi,/is 

whara ia tha apaad of aound. Tha abova aquationa ara baaad on Rafa. (39^ 
and (AO). Becauaa tha anglaa of attack of tha airfoil aactiona ara far 
below the atall angles in all experimental cases p -esantad here, tha abova 
simple relations are applicable. 

The results are shown in Figs. 28 through 37. A comparison of the pre- 
sent theory with the AMI code and with the first experiments^* when using 
the Kocurek and Tangular wake model is presented in Fig. 28. The experi- 
mental result was obtained using a laser velocimeter to measure the bound 
circulation. Although there are differences, especially around x - 0.8, 
the result obtained by using the present theory shows good correlation with 
that of the AMI code. The experimental data show a higher lift peak at the 
blade tip than either theoretical estimation. Tliis is probably caused by 
the existence of a larger wake deformation than estimated using the generalized 
wake roodel'*^, rather than by modeling errors in the two theories. The 
same experimental result^** was analyzed by using the Landgrebe wake model 
as shown in Fig, 29. In addition to the comparison of the present theory 
with the AMI code, the result obtained by the UTRC code is also presented. 

Tlie three theoretical results are very similar. Comparing Fig. 29 with 
Fig. 28, the Landgrebe wake model gives better estimation than the Kocurek 
and Tangular wake model in this experimental case. The difference of lift 
distributions between the two wake models is, however, very small. 
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In addition to tha raaulta ahown in Pig. 29, tha raaulta using tha ! 

pravious rigid waka nodal of tha Local Moaantum Thaory and using blada ala- 
mant thaory with a uniform inducad valocity are shown in Fig. 30. The I 

uniform iducad valocity is given as where tha thrust coefficient 

is equal to that of tha AMI coda. For convanlanca, the lift curve slope t 

i 

j 

of the airfoil section is again given as aquations (72), (73) and (74) in i 

i 

the rigid wake model of the Local Momentum Theory and in blade element 
theory. It is apparent that the present extension of the Local Monumtum 
Tliuory improves the capability of estlractlng the lift distribution near the 
blade tip. 

An analysis of the second experimental case using the Kocurek and 
Tangular wake model is shown in Fig. 31. The experimental data were measured < 

by a flight test using H-34 rotor (Table 4 in Ref. 36). Although the cor- 
relation between the present theory and the AMI code is not as good as tlie 
first experimental case shown in Figs. 28 and 29, it is reasonable. Compared 
with either theoretical prediction, the experimental data show a lower value 
of the lilt distribution at the inner part of the blade. 

Results of another analysis of the second experiemental case^*' using 
the Landgrebe wake model are shown in Fig. 32. The previous results, shown 
in Figs. 28 through 31, were obtained by using the collective pitch input mode 
for all codes, in which the measured value of the collective pitch is given 
to the computer program. In this experimental case, the collective pitch 
input mode showed poor correlation between the UTRC code and the AMI code. 

Therefore, the thrust coefficient input mode, in which the value of the 
thrust coefficient is given to the computer program, was also tried. The 

results using the collective pitch input mode and the thrust coefficient 
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Input Bod* art thown by tht doublt chain lint and by tha broktn lint, 
rtaptotlvtly. Tht thniat ootfflcltnt Input aodt flvta btttar oorrtlatlon, 
but a Btaaux^ablt dlfftrtnot bttwttn tht two ooaputtr oodta ttlll txltta ntar 
tha bladt tip. Btcauaa tht aaat thruat ootfflcltnt la uatd In tht AMI 
code and in tht UTRC coda whtn ualng tht thrmt ootfflcltnt Input aodt, 
the difference between theae two reaulta coata from the difference between 
lifting surface and lifting line, and also cornea from the difference between 
the numerical methods to calculate the thruat coefficient. In idils 
experimental case, the tip vortices generated by the preceding blade are 
closer to the following blade than in the first experimental case shown 
in Figs. 28 through 30. Therefore, tha airloading distribution is more 
sensitive to the tip vortex positions, which depeM on the thrust coefficient. 
This is the primary reason for the poor correlation between the two computer 
codes when using the collective pitch input mode. 

The results of the third experimental case, analyzed by the Kocurek 
and Tangier wake model, are shown in Fig. 33 - The experimental data were 
obtained by a UH-IA helicopter^"^. The correlation between the present 
theory and the AMI code is very good. The correlations between the experiment 
and the present theory or the AMI code are good. It is also apparent 
that the present extension of the Local Momentum Theory shows good 
improvement in the capability of estimating the lift distribution. 

The analyses of the fourth experimental case‘s using the Kocurek and 
Tangier wake model are presented in Figs. 3^ and 35. Fig. 34 shows the 
results using the collective pitch input mode. The experimental data, 


51 


how«v«r, thoKT poor corrolatlon with tb« th*or«tical Mtlmtioiw, Ironically, 
tha pravioua rigid waka modal givaa battar corralation with tha axparimant 
in thia input moda (Rafar to Fig. 16). Whan uaing tha thniat coafficiant 
input moda, tha praaant thaory and tha AMI coda ahow good corralation with 
tha axparimant aa ahown in Fig. 35. Tha airloading diatribution naar tha 
blada tip can ba aatimatad only whan uaing tha diatortad waka modal. 

The fifth axparimantal caaa^® waa analy.ad by tha Kocurak and Tangular 
waka modal. Tha raaulta uaing tha collactiva pitch input moda and tha thruat 
coafficiant input moda ara ahown in Figa. 36 and 37, reapactively. Again, 
only the thruat coefficient input moda of tha praaant thao’*y or of tha AMI 
coda ahow good corralation with tha experiment. Tha correlation between 
tha present thaory and tha AMI coda ara reasonable in both input modes. 

Reviewing analyses of the five experimental cases, correlation between 
the present theory and the AMI code is always good. However, the correlation 
near a blade tip is not as good as that of the inner part of tha blade, nor 
it is as good when the aspect ratio of the blade becomes smaller, or when 
the tip vortices of the preceding blade are closer to the blade under 
consideration. These tendencies of the correlation are due to the simplicity 
of the present model. An Improvement of the attenuation coefficient at the 
upwash region might be effective in obtaining a batter correlation. 

The difference between the results obtained using the AMI code and the 
UTKC code is sometimes greater than the difference between those obtained 
using the AMI code and the present theory. The experimental results fluctu- 
ate around the three theoretical results, and it is difficult to determine 
which computer code gives the best correlation with the experiments. There- 
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for*, th* d*v*lopn«nt of a nor* CMipllcatad nod*l of th* pr*s*nt th*ory, which 
night nak* it poialbl* to obtain a battar corralatlon with th* AMI cod* alona, 
ia not nacAbsary at praaant. 

Th* coaputar tin* r*quir«d for aach of tha thra* codas to obtain tha rasults 
shown Pigs. 29 and 32 is glvan in Tabla 2. Thasa calculations war* conductad 
under the sane condition by using CU07600 at Anas Elasaarch Canter. Ttie 
variation of the computer tin* of th* UTRC cod* in th* case of Fig. 32 is 
due to tha difference of the convergence speeds of the iteration included 
in this code. TTie computer time required by the present theory is less 
than 1/15 of that of the UTRC code, and less than 1/60 of that of the AMI 
code. 


The experimental data of the fourth and the fifth cases are very dif- 
ferent from the theoretical estimation obtained by the present theory and 
by the AMI code when using the collective pitch Input tm>de. This may be 
due to clastic torsion deflection of the blade. When the thrust co- 
efficient input mode is used, however, ; le results obtained from using the 
distorted wake model show the best correlation with the experiments. In 
the other three experimental cases, the results obtained using the present 
theory and the AMI code show reasonable correlation with the experimental 
results, even when the collective pitch input mode is used. 

Conclusion 

The fundamentals of the Local Momentum Theory were presented to make 
clear the boundaries of this theory and to assist to make a computer pro- 
gram using this theory. Tlien, an extension of this theory to include the 
effect of the rotor wake contraction in hovering flight was presented. Tliis 
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«xt«nd«d Locbl Momsntum Theory was comparad with a praacribad waka vortax 
thaory. Tha results Indlcatad that tha extended Local Momentum Theory has 
the capability of achieving a level of accuracy similar to that of a pra- 
scribed wake vortex thaory over wide range variations of rotor geometrical 
parameters. It was also shown chat the analytical results obtained using 
either theory were In reasonable agreement with experimental data. The 
computer time required by the present theory was less than 1/15 of Chat 
of the prescribed wake vortex theory. Consequently, this extended Local 
Momentum Theory is especially effective when calculations are required for 
numerous cases. 
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Appendix A-1 


Recurrence Equation to Give the Induced Velocity Distribution In the 
Case of the Fixed Wing 

In this section, a sanq>le procedure is presented to derive a recurrence 
equation from equation (28). This recurrence equation Is used to solve the 
Induced velocity distribution along the fixed wing span. 


For simplicity, is assumed that the wing has a rectangular plarform, 
no twist angle and the same airfoil section along the span, that is, the 
chord, the geometrical angle of attack and the slope of the lift coefficient 
are constant along the wing sptn. Furthermore, the assumed elliptical wings 
are onesldedly arranged as shown in Fig. 2. In addition to this, the dif- 
ference of the span length is constant between neighboring elliptical wings, 
that is, 'Yj. is constant. 


Substituting j = 3 into equation (28), equation (28) becomes 


* Z 

Also, in the case of j = 1, equation (21) becomes 

= ]',]’■ ft-, 


LA-1-1) 


(.A-1-Z) 


Substituting equation (A. 1-2) into (A.1-f) and performing the integration, 
the following equation can be obtained: 


' t J 


CA-1-3) 
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where 


(A.1-4-) 


In equation (A. 1-3), g(x, y) comes from the Integration and it is generally 
defined by 

‘jU, + S>» X - -Sfri'j) 

(.A-1-S) 


Therefore, g(l,^j^) is given as 


f T(-f- - s>«V^ CA-i-6) 


Similar to the case of J ■ 1, the following can be obtained from 
equations (21) and (28) in the case of j 2: 


AU- = /ev'- 

» ( ipf Jr/ 


hi. ■&(%,%„, bjj 
H>j (A'l-t) 


where 

. f-’i.-, 

"t 

Performing the integration in the above equation and using equation 
(A. 1-4), equation (A. 1-7) becomes 
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wh«r« 


(A-1-10) 


\ 

i 

V.J 

■^fi “ 7*c<) (A-1-io 

As the one-sided arrangement Is assumed here^y^j^ls given by 

7^=a-o-^7/2 (M-ii) 

In summary, the initial value is given by equations (A. 1-3), 

(A. 1-4) and (A. 1-5). The remaining (j i 2) is given by equations 
(A. 1-4), (A. 1-5) and from (A. 1-7) through (A. 1-12). After the values of 
Al^(j ■ l'^'*'*') have been obtained, the induced velocity distribution and 
the lift distribution can be obtained by using equations (27) and (17), 
respectively. 


Appmdlx A-2 


Momentum of g «ch Ellipt ical Wing Qptratlng In « ShMfgd Flow 

In this section, a momentum balance of each elliptical wing operating 
In the sheared flow shown in Fig. 5^. Is discussed. As defined in the sub- 
section of the rotary wing, aselliptical wing means a wing which has an 
elliptical circulation distribution along the span. The horlsontal velocity 
for the elliptical wing is linearly increased from the root to the tip 
along the wing span as given by equation (29). The trailing vortices shed 
from each elliptical wing are straight, perpendicular to the wing span, 
and they are horizontally extended to infinity. The elliptical circulation 
distribution of the wing can be expressed by the following equation using 
a constant Ct : 

nx)= Cr 

= CrV 1- (M-n 

where the relation between^ and x is given as equation (35). Using the 
Blot-Savart relation, the induced velocity inside the wing span is given as 





4-TtR 



dr lx') 
x-x‘ 


(A.Z-J!) 


Substituting equation (A. 2-1) into (A. 2-2) and performing the integration, 
equation (A. 2-2) becomes 
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Substituting Cy * 2 R(|-X^)aI^ from ths abovs aquation into 

(A. 2-1), equation (A. 2-1) is rewritten as 


r(x) = 2R(l-Xi.)^utV/-i* 

The spanwise lift distribution of each elliptical wing, i.e., (/i. «;)«/» 
is given as 

(Uu>)^i^= fUV {A-2-S-) 

Substituting equations (29) and (A. 2-4) into (A. 2-5), equation (A. 2-5) becomes 

( 2 fV,.c R 

x/TTP" 


The total lift of each elliptical wing can be given as 

U* r fV(x)r(x) Rdx (A-Z-l) 

Substituting equations (29) and (A. 2-4) into the above equation, can be 
expressed by 


Li, = 2 m-i, ^ 


(A.2-?; 


1*1;, =/’7i:|f?(i-Xi,y2)^l4,t 

Vi,c~ l^5»«V^ + XiR(i + 


where 
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(A-2-9) 

(A-2-IO) 


Appendix A-3 

R«cncr»nci Equ«t:lon_tP aivg th» Induct _Vftlaclty Pli£rlbatiQn_liL__tkt_£iaA 
of the Rotary Wing 

In this section, a sample procedure Is presented to derive a recurrence 
equation from equations (37) and (39)-(47). The recurrence equation Is used 
to solve the spanwlse Induced velocity distribution of the rotary wing. 

For simplicity, It Is assumed that a rotor blade has a rectangular 
planform, no twist angle and the same airfoil section along the span, that 
is, the chord, the pitch angle and the slope of the lift coefficient are 
constant along the blade span. Furthermore, it Is assumed that the difference 
of the span length is constant between neighboring elliptical wings, that 
is, AX-** Xi/ia constant. 

Using equation (35) and replacing x In equation (47) withj , equation 
(47) becomes 

•'f. .* 




where 


AX * ~~ x.(/ 


(A3-3) 

(A3 -4) 

(4.3- i-; 


and where Is given as equation (33). Performing the Integration in 

equation (A. 3-1), 9 ^ can be given as 




uO 


In the above equation, H(x, y) cornea from the Integration and it is defined 
by 


■* \/l,c 


where g(x, y) is given as equation (A.i-5). Replacing in equation (39) 
with (A. 3-b), substituting i’ » 1 into equation (39) and using equations 
(37), (45) and (A. 3-5) ,AVt can be obtained. 


5 = (Bill - - l/iti. 



a^CAXUi 



CA^}-8) 


where (J, ^ 9^ J,, I are given by equations 

(40), (46), (A. 3-7), (A. 3-3) and (A. 3-4), respectively. 


Similar to i’ ■ 1, the following can be obtained in the case of i' ^ 2: 

(eUa - [1 + H 


where U|», Vjj, H(x, y), and are given by equations (40), 

(46), (A. 3-7), (A, 3-3) and (A. 3-4), respectively. 
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In sunmary, whan known, the Initial value le given by equa- 

tion (A. 3-8), and the remaining ^ X^« (t/|[2) Is given by equation (A. 3-9). 
After the values have been determined, the Induced velocity 

distribution and the lift distribution can be obtained by uslnv equations 
(37) and (38) respectively. The upwash distribution outside the blade tip 
can be obtained by using equation (52). 
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A 

D 

C 

D 

E 

Blade section 

0012 

0012 

0012 

0015 

0012 

Blade twist, 
(deg.) 

-10.902 

-8.0 

-12.0 

0 

0 

Collective pitch Ingle at x 

- 0.75, 
9.41 

8.97 

8.0 

5.3 

Rotor radius, 
R(m) 

1.045 

8.53 

6.67 

0.762 

2.32 

Number of blades, 
b 2 

4 

2 

2 

2 

Solidity, 

0.0464 

0.0621 

0.0369 

0.0637 

0.0974 

Aspect ratio of 

/(R 

a blade, 
13.7 

20.5 

17,3 

10.0 

6.54 

Flapping hxage 

offset, 

0 

0.0357 

0 

0 

0 

Mach number of 

blade tip, 
0.226 

0.583 

0.639 

0.188 

0.360 

Source 

NASA TM 
78615(34) 

NASA TM 
X-952(36) 

TCREC- 

TR-62- 

42(37) 

NACA TN 
2953(26) 

NACA TN 
3688(38) 


Model 

H-34 

Flight 

Test 

UH-IA 

Flight 

Test 

Model 

Model 


Table 1 Rotor parameters and operating conditions 
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CoBq>uter time («ec) 


Code/ Case 

Fig. 29 

Fig. 32 

Present theory 

0.34 

0.44 

AMI code 

23.18 

43.03 

UTRC code 

6.37 

7.02-40.52 


Table 2 Comparison of coirput-T time 



Fig. 1 Lift and induced velocity diatributiona of an elliptical wing'. 





Fig, 2 Decomposition of a fixed wing to hypothetical elliptical wlngsl 
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LIFT FROM BLADE ELEMENT THEORY 
^ LIFT FROM MOMENTUM THEORY 



LIFT DISTRIBUTION 


(j+DTH ELLIPTICAL WING 
jTH ELLIPTICAL WING 
(j-l)TH ELLIPTICAL WING 



T? » y/(b/2) 


Fig. 3 Lift acting on a local segment of a fixed wing. 
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LOCAL MOMENTUM THEORY (ONE-SIDED ARRANGEMENT) 
LOCAL MOMENTUM THEORY (SYMMETRIC ARRANGEMENT) 
LIFTING LINE THEORY BY MULTHOPP20 




(a) LIFTD STRIBUTION 


(b) INDUCED VELOCITY DISTRIBUTION 

Ffj’. 4 I-it t and induced velocity d irft ribut ioiiH l or a rectanju>5 
(AR = 50). 




AIRLOADING DISTRIBUTION 
IS NOT ELLIPTICAL BUT 



Fig. 5 n*compoaltion of a rotor blade to hypothetical elliptical wings. 



Summation of induced velocity distributions of component elliptical 














V 



Fig. 9 Effect of upwash of a preceding blade on the following blade^. 
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Fig. 10 Global flow chart. 
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Fig. 1 
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It 

lU 
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0 



.7 ,8 .9 1.0 


ATTENUATION COEFFICIENT, C 


1 Tlirust coefficient versus uniform attenuation coefficient in hover^. 
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SPANWISE AIRLOADING, i (kg/m) 


PRESENT THEORY 

PRESENT THEORY WITH UNIFORM C* - 0.80 

VORTEX THEORY (ICHIKAWA26) 

A EXPERIMENTAL RESULTS (NACA TN 295326) 



Fig, 16 Spanwlse airloading of a rotor blade in hover^. (p “ 0, b » 2» 
6t * 0 deg.) . 
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SPANWISE AIRLOADING, 8 SPANWISE AIRLOADING, £ 


^ - 0° PRESENT THEORY 

EXPERIMENT (NASA TN D-ISST^^) 


E 

200 




zii A 


X *r/R 

i// = 180° PRESENT THEORY 
= 186° EXPERIMENT (NASA TN D-163?27) 



0 .2 
(b) ^ = 180° 


.4 .6 

x = r/R 


Fig. 17 Spanwise airloading of a rotor blade in forward flight^ (p 
b * 4, &t~ “8.0 deg.). 


A x-0.75 EXPERIMENT (NASA TN 0-1637)27 

x-0.74 PRESENT THEORY 

X - 0.75 RIQIO WAKE VORTEX THEORY^O 

X - 0.75 OISTORTEO WAKE VORTEX THE0RY28 



Fig. 18(a) Azimuth-wise variation of air loading in forward flight^ (w - 0.18, 
b - A, 6|^— 8.0 deg.). 
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TRACE OF INTERSECTION OF BLADE 

UNDER CONSIDERATION WITH TIP 
VORTEX OF PRECEDING BLADE 

Fig. 20 Distribution of angle of attack (deg.)**^ (y » 0.18, b *• 4, 0t=-8.O 



AERODYNAMIC COMPONENT OF Cj 


PRESENT THEORY 

PRESENT THEORY WITH ASSUMED Vq 

VORTEX THEORY** 


i 


1 


j 



2 


cS“ 


0 360 720 1080 

\l/, deg 

(a) THRUST COEFFICIENT 


PRESENT THEORY 

CLASSIC MOMENTUM THEORY'*® 



0 360 720 1080 


i/.deg 

(b) INFLOW RATIO AT r/R « 0.75 

Fig. 22 Comparison of the results obtained by present theory and classical 
momentum theory^ (u ■ 0, b “ 3, 0t “ 0* 0 deg.). 
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LATERAL FLAPPING. 


0 



P 


— — MOMENTUM THEORY INCLUDING A ROTOR 

APPARENT MASS2 

— •*— MOMENTUM THEORY WITH NONUNIFORM 

INDUCED VELOCITY DISTRIBUTION^ 

— -MOMENTUM THEORY WITH UNIFORM 

INDUCED VELOCITY DISTRIBUTI0N32 

Fig. 23 Inclination of tip path plane following a step cyclic pitch input^ 
(u -0,b*4,(?^“0 deg.). 
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LONGITUDINAL FLAPPING, 



, 1 . 8 " 

01 , — L - 4 . 8 “ -iy 

(a) TIME RESPONSE 



(b) INCLINATION OF TIP PATH PLANE 


Fig. 24 Time response to a step cyclic pitch input^ (y » O.IH, h 
-B.O deg.). 
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NON-DIMENSIONAL AIRLOADING 

e/p«2R3 



Fig. 28 Comparison of measured and calculated snanwise lift distributions 
in hover, with Kocurek and Tangular wake model, (b = 2,6t“ -10.9 
deg. ,/^R - 13.7,6'p^^ 9.8 deg.). 
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— PRESENT THEORY (LANDGREBE WAKE) 

— • — PRESCRIBED WAKE THEORY (LANDGREBE WAKE) 

AMI CODE 

— — — PRESCRIBED WAKE THEORY (LANDGREBE WAKE) 

UTRC CODE 

A EXPERIMEN r (NASA TM 78615)3^ 



Fig. 29 Comparison of measured and calculated spanwlse lift distributions 
In hover, with Landgrebe wake model, ( b 2.6t “ -10.9 deg., 

/IR- 13.7, e«.7S- 9.8 deg.). 
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NON-DIMENSIONAL 


PRESENT THEORY (LANDGREBE WAKE) 

PRESENT THEORY (RIGID WAKE) 


BLADE_^I.EIV1ENT THEORY WITH UNIFORM DOWNWASH 

(v * C-|-/2) 

A EXPERIMENT (NASA TM 78615)34 



Fig. 30 ComparlHon of distorted vrake model with rigid wake model or blade 
element theory, (b “ 2,6t ^ -10.9 deg.,/4R>* 13.7, 0o,75* 9.8 deg.). 
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PRESENT THEORY (K-T WAKE) 

PRESCRIBED WAKE THEORY (K-T WAKE) AMI CODE 

A EXPERIMENT (NASA TM X-962)3® 



Fig. 31 Comparison of measured and calculated spanwise lift distributions 
in hover, with Ko'urek and Tangular wake model, (b « 4,6t* -B.O 
deg.,AR * 20. 5,0, vis* 9 .41 deg. ) . 
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PRESENT THEORY (LANDQREBE WAKE); tfo.75 * 9Ar 

PRESCRIBED WAKE THEORY (LANDGREBE WAKE) 

AMI CODE: 00.76 “9.41“ 

— • — PRESCRIBED WAKE THEORY (LANDQREBE WAKE) 
UTRC CODE: 00.75 “ 9.41“ 

. . . PRESCRIBED WAKE THEORY (LANDGREBE WAKE) 
UTRC CODE: 0o.7B“1O>4“ 

A EXPERIMENT (NASA TM X-952)36: 0o.75 * 9.41“ 



Fig. 32 


Comparison of measured and calculated spanwise lift distributions 
in hoverj with Landgrebe wake model , (b - 4,6t “ ~d.O deg., 


(IR - 20.5). 
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— PRESENT THEORY (K T WAKE) 

— • — PRESCRIBED WAKE THEORY (K T WAKE) AMI CODE 

— BLADE ELEMENT THEORY WITH UNIFORM DOWNWASH 

— — - PRESENT THEORY (RIGID WAKE) 

A EXPERIMENT (TCREC-TR-62-42)^^ 



Fig. 33 Comparison of measured and calculated spanwlse lift distributions 
in hover, with Kocurek and Tangular wake model, (b » 2, Ot * ■‘12.0 
deg.,K^“ 17.3, 8.97 deg.). 
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F. 


PRESENT THEORY (K-TWAKE) 

— . — PRESCRIBED WAKE THEORY (K T WAKE) AMI CODE 
A EXPERIMENT (NACATN 2953)26 



r/R 


Fig. 34 Comparison of measured and calculated spanwlse lift distributions 
in hover, collective pitch input mode, (b ■ 2,0t- 0.0 deg. , 

AR- 10.0 deg, 8.0 deg.). 
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— PRESENT THEORY (K T WAKE) ! OqJB ** «•««’ 

»•— PRESCRIBED WAKE THEORY (K T WAKE) AMI CODE 

% JB ■ 8-5S* 

— - - PRESENT THEORY (RlQir- WAKE) ■* Oq.75 * <9 0“ 

A EXPERIMENT (NACA TN 2953)26 ; dojB ” 8,0^ 



r/R 


Fig. 35 


Con^arison of measured and calculated spanwise lift distributions 
in hover, thrust coefficient input mode, (b » 2, et- 0.0 deg. , 


AR - 10.0). 
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PRESENT THEORY (K T WAKE) 

. — prescribed wake theory (K T WAKE) AMI CODE 



r/R 


Fig. J6 Corapariaon of measured and calculated spanwiae lift distributions 
in hover, collective pitch input mode, (b * 2,0i ■ 0.0 deg., 

« b.54, 0O.75* 5.3 deg.). 
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NON-DIMENSIONAL 

AIRLOADING, 


— — PRESENT THEORY (K T WAKE): So.76 “ ^ ^ 
- PRESCRIBED WAKE THEORY (K T WAKE) 
AMI CODE:So .75 “ 4.04* ^ 

A EXPERIMENT (NACA TN 3688) : ,75 - 5.3^ 



0 .5 1.0 


r/R 

Fig. 37 Comparison o£ measured and calculated spanwlse lift distributions 
In hover, thrust coefficient Input mode, (b “ 0.0 deg., 

6.54). 
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